Replica and Glasses 



Giorgio Parisi 

Dipartimento di Fisica and INFN, Universita di Roma La Sapienza 
P. A. Moro 2, 00185 Roma (Italy) giorgio.parisi@romal.infn.it 

February 1, 2008 

Abstract 

In these two lectures I review our theoretical understanding of spin glasses paying a 
particular attention to the basic physical ideas. We introduce the replica method and 
we describe its probabilistic consequences (we stress the recently discovered importance 
of stochastic stability). We show that the replica method is not restricted to systems 
with quenched disorder. We present the consequences on the dynamics of the system 
when it slows approaches equilibrium are presented: they are confirmed by large scale 
simulations, while we are still awaiting for a direct experimental verification. 

1 Introduction 

Many progresses have been done in the last years in the study of glass transitions using the 
replica approach. The basic ideas are those of Kauzmann, Adams, Di Marzio and Gibbs, 
however the use of modern theoretical techniques gives a much better insight. 

The basic hypothesis in this approach is that at low temperature the energy (or better 
the free energy) landscape contains an exponential number of minima. The system may lives 
in one of this mimima, but the entropy of the system will get contribution also from the so 
called configurational entropy or complexity. 

I will concentrate on three points in my lectures. 

• The study of soluble models in which the this scenario is exact. The advantage of 
studying these model are the following: 

— You are able to obtain a precise formulation and define in a clear an precise 
manner the various concepts involved. 

— The resulting picture is much more complex that what you can be naively expect 
and a much better inside is obtained, 

— You also obtain new qualitative predictions. 

— It is possible to positioning the mode coupling computations inside this scenario. 
This is crucial in understand the utility and the limits of mode coupling approach. 

• Model independent predictions for the generalization of the fluctuation dissipation 
theorem in off-equilibrium aging dynamics. The fluctuation dissipation theorem does 
not hold when the system is reaching equilibrium and new relation are present, which 
are model independent: they have been tested in numerical simulation and we are 
eagerly waiting for an experimental test. 



• The construction of first principle analytic computations which allows the computation 
of the property of glass forming quantity at all temperatures (i.e. above and below the 
glass transition) by writing down and solving the appropriate integral equation for the 
correlation function. 

The tools used are mostly the replica theory and generalized two times mode coupling 
equations for the dynamics. 

In these lecture I will not present the detailed proof of the various result: I am trying to 
distill the conclusions of about one hundred papers. I will mostly state the main results, in 
some cases sketch the proof and refer to the original literature for more details. 

2 The basic scenario 

In the nutshell many of the ideas I am going to present are not new: they are already in the 
original papers of Gibbs and Di Marzio. However the comparison of glasses and generalized 
spin glasses, introduced in ref. m allow us to put these ideas in a sharper form and to 
test them in numerical (and eventually real) experiments. In this talk I will not discuss 
the theoretical basis under which this scenario has been derived (i.e. the mathematical tool 
needed to derive the results stated for the generalized spin glasses) but I will concentrate 
the attention of the physical picture. 

The basic ideas are quite simple ^ ^, ^, ^, 0, |^. Let us consider a system of particles 
and let us assume that we can introduce a free energy functional F[p\ which depends on the 
density p{x) and on the temperature. We suppose that at sufficiently low temperature this 
functional has many minima (i.e. the number of minima goes to infinity with the number 
(A^) of particles). Exactly at zero temperature these minima coincide with the mimima of 
the potential energy as function of the coordinates of the particles. Let us label then by 
an index a. To each of them we can associate a free energy and a free energy density 

= F^/N. 

In this low temperature region we suppose that the total free energy of the system can 
be well approximated by the sum of the contributions to the free energy of each particular 
minimum: 

Z = exp{-(3Nfs) = ^exp(-/5iV/„). (1) 

a 

When the number of minima is very high, it is convenient to introduce the function 
A/'(/, T, A^) which is the density of minima whose free energy is equal to /. With this 
notation we can write the previous formula as 

Z = I dfexp{-pNf)Af{f,T,N). (2) 

In the region where is exponentially large we can write 

Ar(/,T,Ar)^exp(ArS(/,r)), (3) 

where the function S is called the complexity or the configurational entropy (it is the contri- 
bution to the entropy coming from the existence of an exponentially large number of locally 
stable configurations). 

The relation (^ is valid in the region fm{T) < f < fuiT). The minimum possible value 
of the free energy is given by fm(T). Outside this region we have that Af{f,T) = 0. It all 
cases known T,{fm(T),T) = 0, and the function S is continuous at fm- 



For large values of we can write 



expi-NPfs) ^ dfexpi-N{Pf - S(/,T)). (4) 

fm. 

We can thus use the saddle point method and approximate the integral with the integrand 
evaluated at its maximum. We find that 

/?/s = mm$(/)^/5r-E(r,T), (5) 

where 

$(/)=/3/-S(/,T). (6) 

(This formula is quite similar to the well known homologous formula for the free energy 
,i.e. Pf = miiiEilSE — S{E)), where S{E) is the entropy density as function of the energy 
density) 

If we call /* the value of / which minimize $(/)• we have two possibilities: 

• The minimum is inside the interval and it can be found as solution of the equation 
P = dH/df. In this case we have 

= pf* - S*, E* = S(/*, T). (7) 

The system may stay in one of the many possible minima. The number of minima at 
which is convenient for the system to stay is exp(A^S*) . The entropy of the system is 
thus the sum of the entropy of a typical minimum and of S*, which is the contribution to 
the entropy coming from the exponential large number of microscopical configurations. 

• The minimum is at the extreme value of the range of variability of /. We have that 
/* = fm, and $ = fm- In this case the contribution of the complexity to the free energy 
is zero. The different states who contribute the the free energy have a difference in 
free energy density which is of order (a difference in total free energy of order 1). 
Sometimes we indicate the fact that the free energy is dominated by a few different 
minima by say the the replica symmetry is spontaneously broken p|, ITO|. 



Form this point of view the behaviour of the system will crucially depend on the free 
energy landscape |lT|, i.d. the function S(/, T), the distance among the minima, the height 



of the barriers among them... Although our final task should be to compute the properties of 
the free energy landscape from the microscopic form of the Hamiltonian, we can tentatively 
assume that the landscape of fragile glasses is similar to that of some soluble long range 
models in presence of quenched disorder |TT . 



I cannot discuss here in details the rationale for this hypothesis; it is also clear that 
it cannot be exact and some differences should be present among the predictions of the 
mean field approximation and the real world. Here I will present the scenario coming from 
mean field, stressing some of the predictions that should have a wider range of validity and 
comparing them with numerical simulations for fragile glasses. 

We can distinguish a few temperature regions. 

• For T > Tf the only minimum of the free energy functional is given by a constant 
density. The system is obviously in the fiuid phase. 
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For Tf > T > To there is an exponentially large number of minima with a non- 



constant density p{x) p, |T2|. There are values of the free energy density such that 
the complexity S is different from zero, however the contribution coming from these 
minima is higher that the one coming from the liquid solution with constant p{x). 

• The most interesting situation happens in the region where To > T > Tk- In this 
region the free energy is still given the fluid solution with constant p and at the same 
time the free energy is also given by the sum over the non trivial minima 0, |^. 

Although the thermodynamics is still given by the usual expressions of the liquid phase 
and final free energy is analytic at To, below this temperature the liquid phase corre- 
spond to a system which at each given moment may stay in one of the exponentially 
large number of minima. It is extremely surprising that the free energy of the liquid 
can be written in this region as the sum of the contribution of the minima, according 
to formula (^). 

The time to jump from one minimum to an other minimum is quite large: it is an 
activated process which is controlled by the height of the barriers which separate the 
different minima. The correlation time will become very large below and for this 
region Td is called the dynamical transition point. The correlation time (which should 
be proportional to the viscosity) should diverge at Tk- The precise form of the this 
divergence is not well understood. It is natural to suppose that we should get divergence 
of the form exY>{^/ {T — TkY) for an appropriate value of v [|T^|, whose reliable analytic 
computation is lacking |T^ . 

The equilibrium complexity is different from zero (and it is a number of order 1) 
when the temperature is equal to Tq and it decreases when the temperature decreases 
and it vanishes linearly at T = T^- At this temperature (the so called Kauzmann 
temperature) the entropy of a single minimum becomes to the total entropy and the 
contribution of the complexity to the total entropy vanishes. 

• In the region where T < Tk the free energy is dominated by the contribution of a few 
minima of the free energy having the lowest possible value. Here the free energy is no 
more the analytic continuation of the free energy in the fluid phase. A phase transition 
is present at Tk and the specific heat is discontinuous here. 

The free energy landscape is quit usual; we can try do present the following pictorial 



interpretation, which is a rough simplification |T5[. At temperature higher than the 
system stays in a region of phase space which is quite fiat and correspond of a minimum of 
the total free energy. On the contrary below To the phase space is similar to the one shown 
pictorial in fig. 1. The region of maxima and minima is separated by the region without 
barriers by a large nearly flat region. The minima in the region at the left are still present 
also when Tf > T > T^,, but they do not correspond to a global minimum. 

At temperatures higher than the system at thermal equilibrium states in the right 
region. When the temperature reaches To the system arrives in the fiat region. Here the 
potential is fiat and this causes a more or less conventional Van Hove critical slowing down 
which is well described by the well known mode coupling theory [|16| (which is exact in the 



mean field approximation). The mode coupling theory describes the critical slowing down 



which happens near To |T7| 



In the mean field approximation the height of the barriers separating the different minima 
is infinite and the temperature To is sharply defined as the point where the correlation 
time diverge. In the real world activated process (which are neglected in the mean field 
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Figure 1: The qualitative dependence of the free energy as function of the configuration 
space in the region relevant for the dynamical transition, i.e. for T < T^. 

approximation and consequently in the mode coupling theory) have the effect of producing 
a finite (but large) correlation time also at To. The precise meaning of the dynamical 
temperature beyond mean field approximation is discussed in details in |^ 

When the temperature is smaller that we must be more precise in describing the 
dynamics of the system. Let us start from a very large system (of particles) at high 
temperature and let us gradually cool it. We find that it should go at equilibrium in the 
region with many minima. However coming from high free energy (from the right) it cannot 
enter in the region where are many maxima and minima if we wait a finite amount of time 
(the time to crosses the barriers diverges as exp{AN). If we do not wait an exponentially 
large amount of time the system remains confined in the fiat region. In this case the so 
called dynamical energy, 

Ed = lim lim E(t,N), (8) 

is higher that the equilibrium free energy. The situation is described in fig. 2. 

Below Td the system is trapped in metastable states when cooled. The time needed 
to escape from these states diverges when goes to infinity in the mean field approach 
where activated processes are forbidden. Of course the difference of the static and dynamic 
energy is an artifact of the mean field approximation if we take literarily the limit t ^ oo in 
the previous equation because as matter of fact there are no metastable states with strictly 
infinite mean life. However it correctly describe the situation on laboratory times, where 
metastable states are observed. 

In real systems, beyond the mean field approximation, the height of the barriers is finite 
also below Td and the mean life of the metastable states is finite, albeit very large. Some 
mechanisms have been described who imply a divergence of the correlation time in real 
systems at the Kauzmann temperature [||, ^, |T^. The conventional glass temperature, i.e. 
the temperature at which the microscopic correlation time becomes macroscopic (e.g. of 
order of the minute) is between the two temperatures (i.e. Td > Tq > T^)- 

It should be clear that in this framework the dynamical temperature is not so well 
defined and it correspond to a crossover region below it the dynamics is dominated by 
activated processes 

In the mean field approximation there are very interesting phenomena that happen below 
To when the system is cooled from the high temperature phase. These phenomena are related 
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Figure 2: The qualitative behaviour of the equihbrium and of the dynamical energy as 
function of the temperature. 



to the fact that the system does not really go to an equilibrium configuration but wanders 
in the phase space never reaching equilibrium. The phenomena are the following: 

• The energy approaches equilibrium slowly when the system is cooled from an high 
energy configuration. In other words for large times we have 

E{t, T) = Ed{T) + B{T)t-^'^^\ (9) 

where the exponent A(T) does not vanish linearly at zero temperature as it should 
happens for an activated process. 

• Aging is present, i.e. the correlation functions and the response functions in the region 
of large time do depend on the story of the system [0, [1^, . 

• In the region where aging is present the hypothesis at the basis of the fluctuation 



dissipation theorem are no more valid. New generalized relations are satisfied |^ 
which replace the fluctuation dissipation theorem. 

• In the region where the diffusion constant is zero, a new phenomenon, logarithmic 
diffusion, is present. 

All these phenomena are well known also in the case of spinodal decomposition but have a 
more general validity. 



3 The Random Energy Model 
3.1 The definition of the model 

The Random Energy Model is the simplest model for glassy systems. It have various ad- 
vantages: it is rather simple (its properties may be well understood with intuitive arguments, 
which may become fully rigorous) and display very interesting and new phenomena. 

The Random Energy Model is defined as following. There are Ising spins (aj, i = 1, A^) 
which may take values ±1; the total number of configurations is equal to M = 2^ and they 
can be identified by a label k in the interval 1 — M. 



Generally speaking the Hamiltonian of the system is given when we know the values of 
the energies for each of the M configurations of the system. Usually one writes explicit 
expression for the energies as function of the configuration; on the contrary here we assume 
that the values of the are random, with a probability distribution p{E) which is Gaussian: 

p(E)^expi-^). (10) 
The partition function is simply given by 

Ze= exp(-/3Efc)= f p{E)exp{-PE), (11) 

k=l,M 

p{E)^ Y: S{E-E,). 

k=l,M 

The value of the partition function and of the free energy density (/f = — In^Z^)/ (NP)) 
depends on S, i.e. all the values of the energies E^. We would like to prove that when 
N ^ oo the dependance on £ of fs disappears with probability 1. If this happens, the most 
likely value of fs coincides with the average of fs, where the average is done respect to all 
the possible values of the energy extracted with the probability distribution eq. (|r^) . 

The model is enough simple to be studied in great details; exact expressions can be 
derived also for finite N. Here we sketch the results giving a plausibility argument. 



3.2 Equilibrium properties of the model 

The crucial observation is the following. The probability of finding a configuration of energy 
E is 

KiE) ^ 2^exp(-|^) = exp(iV(ln(2) - K')), (12) 

where e = E/N is the energy density. It is reasonable to assume (and it is confirmed by a 
detailed computation) that in the case of a generic system (with probability 1 when —>■ oo) 
no configurations are present in the region where A/'o(-E) << 1, i.e. for 

< = 21n(2). (13) 

We can thus write for the generic choice of the energies (S): 

p{E) ^ Af{E) = No{E)e{El - E^). (14) 

The partition function can be written as 

J Af{E)exp{-(3E). (15) 

Evaluating the integral with the saddle point method one finds that in the high temperature 
region, i.e. 

P<Pc^e;\ (16) 

the internal energy density is simple given by —jS. This behaviour must end somewhere 
because we know that the energy is bounded form below also when (3 —>■ oo. Indeed in 
the low temperature region one finds that the integral is dominated by the boundary region 
E ^ Ec and the energy density is exactly given by — Cc- 

The entropy density is positive in the high temperature region, vanishes at (3c and remains 
zero in the low temperature region. It follows that in the high temperature region an 
exponentially large number of configurations contributes to the partition function, while 
in the low temperature region it is possible that the probability is concentrated on a finite 
number of configurations. 



3.3 Properties of the low temperature phase 

It is worthwhile to study in more details the structure of the configurations which mostly 
contribute to the partition function in the lower temperature phase. At this end it is useful 
to sort the configurations with ascending energy. We rename the configurations and we 
introduce new labels such that {Ek < Ei for k < i). 
It is convenient to introduce the following quantity: 

_ exp(-/jEfc) 

Wk = ^ ■ (17) 

We have obviously that 

Wk = l (18) 

A detailed computation shows |23|, ^, |T0[ that in the low temperature region (i.e. /3 > Pc) 
the previous sum is dominated by the first terms. Indeed 
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Y Wk = l- 0{L-^), A = (19) 

where 



k=l,L ^ 



m = ^. (20) 

In the same region the sum in the following equation is convergent and its average value is 
given by 

E wl = l~m. (21) 

k=l,2^ 

Generally speaking one finds that one can introduce the quantities Fk such that 

Wk cx exp{-(3Fk) (22) 

Here the variables Fk coincide with the total energy (not the energy density!) apart form 
an addictive constant. Their probability distribution at the lower end (which is the relevant 
region for thermodynamics in the low temperature region) can be approximated as 

P{F) ^ expipmF). (23) 

In this model everything is clear: in the high temperature region the number of relevant 
configurations is infinite (as usual) and there is a transition to a low temperature region 
where only few configuration dominates. 

This phenomenon can be seen also in the following way. We introduce a distance among 
two configurations a and 7 as 

d^(a,7)^ ^-^'"^f (24) 
Sometimes it is convenient to introduce also the overlap gdefined as 

g(«,7) = ^^^^^^ = l-ci^(«,7). (25) 

The distance squared is normalized in such a way that it spans the interval — 2. It is equal 

to 
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0, if the two configuration are equal (g = 1). 



• 1, if the configuration are orthogonal (g = 0). 
. 2, ifaf = -a7 (g = -l). 

It is convenient to introduce the function Q{d) and P{q), i.e. the probability that two 
equilibrium configurations are at distance d or overlap q respectively. We find 

• For T > Tc 

Qid) = 5id-1), P{d)=5{q). (26) 

• For T < Tc 

Q{d) = {l-A)5{d~l) + A5{d), P{d) = {l- A)5{d) + A5{q-l). (27) 

where A is equal to X]a:=i,2^ ^fc- The average of A over the different realizations of 
system is equal to 1 — m. 

As soon as we enter in the low temperature region, the probability of finding two equal 
configurations is not zero. The transition is quite strange from the thermodynamic point of 
view. 

• It looks like a second order transition because there is no latent heat. It is characterized 
by a jump in the specific heat (which decreases going toward low temperatures). 

• It looks like a first order transition. There are no divergent susceptibilities coming form 
above of below (which beyond mean field theory should imply no divergent correlation 
length). Moreover the minimum value of d jumps discontinuously (from 1 to 0). 

• If we consider a system composed by two replicas (cr^ and cr^) and we write the 
Hamiltonian 

H{a\ a^) = H{a^) + H{(t^) + Ned'^{ai, a^) (28) 

the thermodynamics is equal to that of the previous model (apart a factor 2) for e = 0, 
but we find a real first order thermodynamic transition, with a discontinuity in the 
internal energy, as soon as e > 0. The case e = is thus the limiting case of real first 
order transitions. 

These strange characteristics can be summarized by saying that the transition is of order 
one and half, because it share some characteristics with both the first order and the second 
order transitions. 

It impressive to note that the thermodynamic behaviour of real glasses near Tc is very 
similar to the order one and half transition of REM. We will sea later that this behaviour is 
typical of the mean field approximation to glassy systems. 

3.4 Dynamical properties of the model 

The dynamical properties of the model can be easily investigated in a qualitative way. In- 
teresting behavior is present in the region where the value of is large with respect to the 
time. Different results will be obtained for different definition of the dynamics if we consider 
some rather artificial form of the dynamics. 

Let us first consider a single spin flip dynamics. In other words we assume that in 
a microscopic time scale scale, which for simplicity we consider of order unit, the system 
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explore all the configurations which differ from the original one by a single spin flip and 
goes to one of them (or remain in the original one) with probability which is proportional 
to exp{—j3H). This behaviour is typical of many dynamical process, like Glauber dynamics, 
monte Carlo, heath bath. 

In this dynamical process each configuration C has N nearby configurations C to explore. 
The energies of the configurations C are uncorrelated to the energy of C, so that they would 
be of order one in most of the case. The lowest energy of the configurations C would be 
of order — (A^ In(A^))^/^, the corresponding energy density (— (ln(A^)/A^)^/^) vanishes in the 
large limit. 

If the configuration C has an energy density e less that zero, but greater than the equi- 
librium energy, the time needed to do a transition to a configuration of lower energy will be, 
with probability one, exponentially large. 

For short times a configuration of energy e will be completely frozen. Only at larger times 
it may jump to a typical configuration of energy zero. At later times different scenarios are 
possible: the configuration comes back to the original configuration of of energy e or, after 
some wandering int the region of configurations of energy density ~ 0, it fells in an other deep 
configuration of energy e'. A computation of the probabilities for these different possibilities 
has not yet been done, although it should not too difficult. 

The conclusions of this analysis are quite simple. 

• If we start from a random configuration, after a time which is finite when A^ — > oo, the 
system goes to a configuration whose energy is of order — In(A^)^/^ and stops there. 

• If we start from a random configuration and we study the system at exponentially 
large times the system will reach an energy density which may be different from the 
original one. 

4 Models with partially correlated energy 
4.1 The definition of the models 

The random energy model (REM) is rather unrealistic in that it predicts that the energy 
is completely upset by a single spin flip. This feature can be eliminated by considering a 
more refined model, the so called p-spins models P^ , p3[ , in which the energies of nearby 
configurations are also nearby. We could say that energy density (as function of the con- 
figurations) is not a continuous function in the REM, while it is continuous in the p-spins 
models, in the topology induced by the distance (eq. ( p^ ). In this new case some of the 
essential properties of the REM are valid, but new features are present. 

The Hamiltonian we consider depends on some control variables J, which have a Gaussian 
distribution and play the same role of the random energies of the REM and by the spin 
variable a. For p = 1, 2, 3 the Hamiltonian is respectively 

H'Aa) = Y: J^^^ (29) 

i=l,N 

Hj{<y) = Ak'^i'^k (30) 

i,k=l,N 

I 

i,k,l=l,N 
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where the primed sum indicates that all the indices are different. The variables J must have 
a variance of 0{N^^~'^^/'-^) in order to have a non trivial thermodynamical limit. 

It is possible to prove by am explicit computation that if we send first N oo and later 
p —>■ oo, one recover the REM. Indeed the energy differences corresponding to one spin flip 
are of order p for large p (they ar order in the REM, so that in the limit p ^ oo the 
energies in nearby configurations become uncorrelated and the REM is recovered. 



4.2 Equilibrium properties of the models 

The main new property of the model is the correlation of energies. This fact implies that 
if C is a typical equilibrium configuration, all the configurations which differ from it by a 
finite number of spin flips will also have a finite energy. The equilibrium configurations are 
no more isolated (as in REM), but they belongs to valleys, such that the entropy restricted 
to a single valley is proportional to and it is and extensive quantity. 

The thermodynamical properties at equilibrium can be computed using the replica me- 
thod [1^ . Let us discuss firstly what happen for p > 2 In the simplest version of this 
method ||2^, ^ one introduces the typical overlap of two configurations inside the same 
valley (sometimes denoted by qEA)- Something must be said about the distribution of the 
valleys. Only those which have minimum free energy are relevant for the thermodynamics. 
One finds that these valleys have zero overlap and have the same distribution of free energy 
as in the REM 

P{F) oc exp(/3m(F - Fq)). (31) 

Indeed the average value of the free energy can be written in a self consistent way as function 
of m and q {f{q,Tn)) and the value of these two parameters can be found as the solution of 
the stationarity equations: 

dm dq ^ ^ 

The quantity q (which would be 1 in the REM) is here of order 1 — exp{—A(3p) for large 
p, while the parameter m has the same dependence of the temperature as in the REM, i.e. 1 
at the critical temperature, and a linear behaviour al low temperature. The only difference 
is that m is no more strictly linear as function of the temperature. 

The thermodynamical properties of the model are the same as is the REM: a discontinuity 
in the specific heat, with no divergent susceptibilities. Let us first recall the description of 
these glassy systems at equilibrium according to the predictions of the replica theory. We 
consider a systems and we denote by C a generic configuration of the system. For simplicity 
we will assume that there are no symmetry in the Hamiltonian, in presence of symmetries 
the arguments must be slightly modified. It is useful to introduce an overlap q{C,C'). There 
are many ways in which an overlap can be defined; for example in spin system we could 
define 

9 = , (33) 

being the total number of spins or particles and a and r are the two spin configurations. 
In a liquid a possibility is given by 

Ei=i,N Ek=i,N f{x{t} - y{k)) 
g = , (34) 

where / is a function which decays in a fast way at large distances and is substantially 
different from zero only at distances smaller that the interatomic distance {x and y are the 
two configurations of the system). 
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In the high temperature phase for very large values of the probability distribution of 
the overlap {P]\f{q)) is given by 

PN{q)^5{q-q*). (35) 

In the low temperature phase PN{q) depends on (and on the quenched disorder, if it is 
present). When we average over we find a function P{q) which is not a simple delta 
function. In all known case one finds that 

P{q) = amS{q - qm) + o-m^Q - 1m) + p{q), (36) 
where the function p{q) does not contains delta function and its support is in the interval 

[qm,qM]- 

The non triviality of the function P{q) (i.e. the fact that P{q) is not a single delta 
function and consequently q is an intensive fluctuating quantity) is related to the existence 
of many different equilibrium states. Moreover the function PN{q) changes with N and its 
statistical properties (i.e the probability of getting a given function Pn{q)) can be analytically 
computed 0, . 

In this equilibrium description a crucial role is given by the function x{q) defined as 

x{q) = r P{q')dq' . (37) 
•Jqm 

In the simplest case the function p{q) is equal to zero, i.e the function P{q) has only two 
delta functions without the smooth part. In this case, which correspond to one step replica 
symmetry breaking, there are many equilibrium states, labeled by a, and the overlaps among 
two generic configurations of the same state and of two different states are respectively qM 
and qm- The probability of finding a state with total free energy / is proportional to 

exp(m/5(/-/H)), (38) 

where fn is a reference free energy and m is the value of x{q) in the interval [qm, qAi]- 

In the more complicated situation where the function p{q) is non zero, couples of different 
states may have different values of the overlaps. The joint probability distribution of the 
states and of the overlaps can be described by formulae similar to eq. (|38|) , but more complex 

i- 

Although these predictions are quite clear, it is not so simple to test them for many 
reasons: 

• They are valid at thermal equilibrium, a condition that is very difficult to reach for 
this kind of systems. 

• Experimentally is extremely difficult to measure the values of the microscopic variables, 
i.e all the spins of the system at a given moment. These measurements can be done 
only in numerical simulations, where the observation time cannot be very large. 

A very important progress has been done when it was discovered that the function X, 
which describes the violations of the fluctuation dissipation theorem, is equal to the function 
X which is relevant for the statics. This equality is very interesting because function X{C) 
can be measured relatively easily in off-equilibrium simulations ||21|| . 

The temperature dependence of the function X{C) (or equivalently x{q)) is interesting 
also because rather different systems can be classified in the same universality class according 
to the behaviour of this function. It has been conjectured long time ago that the equilibrium 
properties of glasses are in the same universality class of some simple generalized spin glass 
models ^ |7|. 
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4.3 The Free energy landscape 

It would be interesting to characterize better the free energy landscape of the model, espe- 
cially in order to understand the dynamics. Indeed we have already seen that in the REM 
the system could be trapped in metastable configurations. Here the situation is more com- 
plicated. Although the word valley has a strong intuitive appeal, we must first define what 
a valley is in a more precise way. 

There are two different (hopefully equivalent) definitions of a valley: 

• A valley is a region of configuration space separated by the rest of the configuration 
space by free energy barriers which diverge when N ^ oo. More precisely the system, 
in order to go outside a valley by moving one spin at once, must cross a region where 
the free energy is higher that of the valley by a factor which goes to infinity with N. 

• A valley is a region of configuration space in which the system remains for time which 
goes to infinity with A^. 

The rationale for assuming that the two definitions are equivalent is the following. We expect 
that for any reasonable dynamics in which the systems evolves in a continuous way (i.e. one 
spin fiip at time), when it goes from a valley to an other valley, the system must cross a 
configuration of higher free energy and therefore the time for escape from a valley is given 

by 

T ~ To exp(/3AF) (39) 

where AF is the free energy barrier. 

It is crucial to realize that in infinite range models there can be valley which have an 
energy density higher that that of equilibrium states. This phenomenon is definitely not 
present in short range models. No metastable states with infinite mean life do exist in 
nature. 

Indeed let us suppose that the system may stay in phase (or valleys) which we denote as 
A and B. If the free energy density of B is higher than that of A, the system can go from 
i? to A in a progressive way, by forming a bubble of radius R of phase A inside phase B. If 
the surface tension among phase A and B is finite, has happens in any short range model, 
for large R the volume term will dominate the free energy difference among the pure phase 
B and phase B with a bubble of A of radius R. This difference is thus negative at large 
R, it maximum will thus be finite. In the nutshell a finite amount of free energy in needed 
in order to form a seed of phase A starting from which the spontaneous formation of phase 
A will start. For example, if we take a mixture of H2 and O2 at room temperature, the 
probability of a spontaneous temperature fiuctuation in a small region of the sample, which 
lead to later ignition and eventually to the explosion of the whole sample, is greater than 
zero (albeit quite a small number), and obviously it does not go to zero when the volume 
goes to infinity. 

We have two possibilities open in positioning this mean field theory prediction of existence 
of real metastable states: 

• We consider the presence of these metastable state with infinite mean life an artefact 
of the mean field approximation and we do not pay attention to them. 

• We notice that in the real systems there are metastable states with very large (e.g. 
much greater than one year) mean life. We consider the infinite time metastable states 
of the mean field approximation as precursors of these Gnite states. We hope (with 
reasons) that the corrections to the mean field approximation will give a finite (but 
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large) mean life to these states (how this can happen will be discussed in the next 
section) . 



Here we suppose that the second possibility is the most interesting and we proceed with 
the study of the system in the mean field approximation. The strategy for investigate the 
properties of these metastable states consists in considering systems with R replicas (two or 
more) of the same system with Hamiltonian given by 



(40) 



r=l,_R 



Different replicas may stay at different temperature. The quantities e are just Legendre 
multipliers needed to enforce specific value of the the overlaps q. In this way (let us consider 
for simplicity the case where all temperature are zero) we find (after a Legendre transform) 
a free energy density as function of the q. 

Let us consider for simplicity the case where we set 



q for r = 2,R 



(41) 



{qr^r is identical equal to 1) and the others q are left free . A simple computation show 
the final free energy density that we obtain (let us call //?((?) is given by 



fniq) - Rf 



lim 



In (E. (E.%(^T,r)-g)) 



(B-l) 



lim 

N^oo 



In < P^(g)(«-i) > 



f3N 



(42) 



where / is the unconstrained free energy. 

A particular case, which is very interesting, is given by the limit R 

N^oo PN OR 
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(43) 



The potential W{q) has usually a minimum at g = 0, where W{0) = 0. It may have a 
secondary minimum at q = qo- We can have three quite different situations 

• W{qD) = 0. This happens in the low temperature region, below T^., where we can put 
two replicas both at overlap and at overlap qEA without paying any prize in free 
energy. In this case qo = qEA- 

• W^qo) > 0. This happens in an intermediate temperature region, above Tc, but below 
To, where we can put one replica a at equilibrium and have the second replica r in 
a valley near it. It happens that the internal energy of both the a configuration (by 
construction) and of the r configuration are equal to the equilibrium one. However the 
number of valley is exponentially large so that the free energy a single valley will be 
smaller. One finds in this way that W{qD) > is given by 



W{q_ 



D 



InTV; 

N 



(44) 



where Me is the average number of the valleys having the equilibrium energy [pi E9 
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• At T > To the potential W{q) has only the minimum at g = 0. The quantity qn 
cannot be defined and no valley with the equilibrium energy are present. This is more 
or less the definition of the dynamical transition temperature Tp. A more careful 
analysis shows that for < T < Ty there are still valleys with energy less than the 
equilibrium one, but these valleys cover a so small region of phase space that they are 
not relevant for equilibrium physics. 

It is also possible to study the properties of the free energy for i? 7^ 1 we can force the a 
configuration to be not an equilibrium one and in this way we control the properties of the 
valleys having an energy different than the equilibrium one. This method can give rather 
detailed information on the free energy landscape, which I do not have time to discuss in 
details and which have not yet yet fully studied. 

The most interesting result is that for T < the entropy of the system can be written 

as 

S = Sv + W (45) 

where Sy is the entropy inside a valley and W is the configurational entropy, or complexity, 
i.e. the term due to the existence of an exponentially large number of states. The W 
contribution vanishes at Tc and becomes exactly equal to zero for T < Tc . 

In the REM limit {p —>■ 00) the temperature goes to infinity. In this limit the third 
region does not exist. Therefore the dynamical transition is a new feature which is not 
present in the REM. 



5 Static and dynamic properties 

The simplest way to study the dynamics of the problem is to consider a system which evolves 
according to Langevin equation of to some sort of Glauber dynamics. For example we can 
suppose that: 

daj SH , . , , 

-7^ = + V^it), (46) 

at ai 

where t] is an appropriate white noise. 

In this model is convenient to introduce the single site correlation function (C) and the 
response (G) function of the times. One finds that they can be defined as 

C{h,t2)=<a,{h)a,{t2)> (47) 

G(t„t,)=<^> (48) 
hi{t2) 

where h{t) is an external magnetic field. 

If the systems is at equilibrium (or in a metastable state), the correlation and the response 
functions will depend only on the time difference. If the system is out of equilibrium these 
functions will depend in a non trivial way from both the arguments. 

In both cases one can write down closed equations for the correlation functions in the 
case where the size N goes to infinity at fixed times . These equations have a rather 
complex structure. We could discuss two different regimes: 

1. We start at time zero from an equilibrium configuration. 

2. We start at time zero from a random configuration. 
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In the first case we find that, if we approach the dynamical temperature from above, the 
correlation time diverges a a power of T — To, and the usual analysis of the mode coupling 
theory can be done in this region. Mode coupling theory is essentially correct in this region. 

In the second case one finds that below the dynamical transition, the energy does not go 
anymore to the equilibrium value, but it goes to an higher value, which in some cases can 
be computed analytically. 

The phenomelogy is rather complex, the aging properties of the systems are particularly 
interesting, but we cannot discuss them for lack of space. It is interesting to note that 
the mode coupling theory become exact in the mean field theory and describes what hap- 
pens nearby the dynamical phase transition at Tq, which is a temperature higher that the 
equilibrium transition temperature at T^. 

Here we would like to understand the dynamical properties of the model starting from 
results which we have already derived on the free energy landscape without having to com- 
pute explicitly the solution of the equation of motion for the correlations C and G. This 
can be done by assuming the this very slow dynamics is an activated process dominated by 
barrier crossing and that the height of the barriers can be computed using the method of 
the previous sections. 

Some of the question we ask are the following: 

• How long does a system remains in the same valley? 

• If the system is out equilibrium at time 0, which is the asymptotic value of the energy 
at large times? 

The results are similar to those obtained by the explicit dynamical analysis of the previous 
section. We obtain also more detailed information on the region where the time scale is 
exponentially large. This region can be studied only with very great difficulties using the 
dynamical equations. 

In the first case (i.e. the system start from a thermalized configuration) below Tp, the 
system is confined to a valley. A first estimate of the time needed to escape from a valley 
can be obtained as follows. We introduce the parameter 

g(t) ^ (49) 

We assume that on the large time scale the evolution of q is more or less the same of 
a system with only one degree of freedom with potential W{q). This assumption is not 
completely correct, but it is likely to be enough to give a first estimate (to be refined later) 
of the escape time. Therefore it is reasonable to assume that in a first approximation the 
the time needed to escape from the valley is given by 

r = Toexp(iV/5AW^), (50) 

where AW^ is the difference in free energy among the minimum at qd and the maximum al 
lower values of q. 

For large time (but not exponentially divergent with A^) we have that 

lim q{t) = qo- (51) 

t— ♦oo 

Eventually the system will escape from a valley for times greater that r. 

In the other situation (the system is originally out equilibrium) things are more compli- 
cated. Metastable valleys of high energies do exist, however it is not clear that a system 
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cooled from an high temperature region must be trapped in the valley with high energy. 
Indeed it will be trapped in the valley with largest attraction domain, but the properties 
of the valley cannot be found if we do not use the dynamics in an explicit way. Some edu- 
cated conjectures can be done, but the question has never been investigated in details. The 
problem is not easy, because the system is strongly out of equilibrium. 

A better understanding comes if we consider the case in which the system is slowly cooled, 
from high to low temperature. We consider here that case of slow, not ultra slow, cooling, i.e. 
the time scale if fixed when N goes to infinity. In this case it is reasonable to assume that the 
system frozen in a valley at the dynamical temperature and we have to follow the energy of 
that valley when we cool. That can be done by using the formulation with different replicas 
at different temperatures. In this way one finds that the system is frozen in a configuration 
which is near to the configuration at the dynamical temperature. 

A detailed comparison of the results for the energy of the metastable states obtained by 
doing different assumption is still lacking, but is should be not too difficult. 



6 Systems without quenched disorder 

Apparently the previous discussions are restricted to systems where quenched disorder is 
present. The requirement of quenched disorder would limit ourselves in the applications of 
the replica method and one would cut all those systems, like glasses, which have a transla- 
tional invariant Hamiltonian and where no quenched disorder is present. 

This prejudice (on the need of quenched disorder( was so strong that it took a few years 
to realize that the replica method can also applied to system without quenched disorder. 

There are many facts that clearly indicate the possibility of applying the replica methods 
to non-random systems. 

• In the infinite range case there are pairs of systems with Hamiltonian respectively Hq 
and H, where Hq contains quenched disorder and no disorder is present in H, such 



the high temperature expansion for the two systems coincide It is natural to 

suppose that the free energies of the two system are identical at all temperatures, so 
that replica symmetry breaking can be applied to both. 

It is possible to look for replica symmetry breaking in the expression for the free 
energy systems without disorder (e.g. soft spheres) inside a given approximation (e.g. 
hypernetted chain) and find out that replica symmetry breaks at low temperature [^]. 

In the replica method we can introduce coupled replica potentials in order to 
characterize the phase space of the system and these potentials can also be computed 
for non-random systems, obtaining the same results as for random systems |]2S|. This 



may be done analytically for soft spheres using the same approximation as before 

• It is now clear that the replica method may be applied any stochastic stable system. 
Indeed stochastic stable systems are the limit of disorder systems where the replica 
method can be apphed without problems. Systems without quenched disorder may 
be stochastically stable if the free energy is computed using the Cesareo limit (i.e. 
averaging over A^). 

This new perspective allows us to use the replica method in systems quite different from 
the usual one, e.g. structural glasses, where no quenched disorder is present. 
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7 Fluctuation and dissipation relations in aging dy- 
namics 



7.1 Theoretical considerations 

When we suddenly decrease the temperature in an Hamihonian system, many interesting 
phenomena happen if the initial and the final temperatures correspond to different phases. 
When the low temperature phase can be characterized by a simple order parameter (e.g. the 
magnetization for ferromagnets) we find the familiar phenomenon of spinodal decomposition 
characterized by growing clusters of different phases. There is a dynamical correlation length 
(i.e. the size of the clusters), which increases as a power of the time t after the quench, and 
the energy approaches equilibrium with power like corrections (e.g. E(t) ~ E^^ + 
In this region aging phenomena are also present 

The situation is more intriguing in the case of structural glasses and spin glasses where 
the low temperature phase cannot be characterized in term of a simple order parameter. 
Remarkable progresses in understanding the off-equilibrium dynamics and its relations to 
the equilibrium properties has been done by noticing that a crucial off-equilibrium feature is 
the presence of deviations from the well known equilibhum fluctuation-dissipation relations. 
On the basis of analytic results for soluble models it has been conjectured that we can deflne 
a function X{C), C being an autocorrelation function at different times |]^, |2D|, [T^. This 
function characterizes the violations of the fluctuation-dissipation theorem (which is correct 
only at equilibrium). It is remarkable that (at least in the case of spin glasses) the function 
X{C) is equal to the function x{q) (g being the overlap of two spin conflgurations) which 
plays a central role in the equilibrium computation of the free energy 0. 

Let us be more precise. We concentrate our attention on a quantity A{t). We suppose 
that the system starts at time t = from an initial condition and subsequently it remains at 
a flxed temperature T. If the initial conflguration is at equilibrium at a temperature T' > T, 
we observe an off-equilibrium behaviour. We can deflne a correlation function 



and the response function 



C{t,t^) = {A{t^)Ait + t^)) (52) 
6{A{t + t^)) 



G{t, tui) = 



(53) 

e=0 



where we are considering the evolution in presence of a time dependent Hamiltonian in which 
we have added the term / dte(t)A{t). 

The usual equilibrium fluctuation-dissipation theorem (FDT) tells us that 



where 



G^\t) = -P^^, (54) 

G^\t)= lim G{t,t^), C^\t)= lim C{t,tyj). (55) 
It is convenient to deflne the integrated response: 

R{t,t^)= dTG{t-r,t^ + T), R^\t)= lim R{t,t^), (56) 

Jo tyj—*(Xj 

R(t, tyj) is the response of the system at time t + t^o to a fleld acting for a time t starting at 
tw The usual FDT relation becomes 

R^i^t) = (3{C^'^{t) - C"='?(0)). (57) 
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Figure 3: The response R times T versus C at T = 0.7 for the three dimensional Ising spin 
glass [10]. The curve is the prediction for function R{C) obtained from the equilibrium data. 
The straight line is the FDT prediction. We have plotted the data of the two runs: = 10^, 
and t,„ = 10^. 



The off-equilibrium fluctuation-dissipation relation states that the response function 
and the correlation function satisfy the following relation for large t^,: 



R{t,t^) ^ (3 / X{C)dC. 



(58) 



If we plot R versus (3C for large tyj the data collapse on the same universal curve and the 
slope of that curve is —X{C). The function X{C) is system dependent and its form tells us 
interesting information, in the case of three dimensional spin glasses. 
We must distinguish two regions: 

• A short time region where X{C) = 1 (the so called FDT region) and C belongs to the 
interval I (i.e. Ci < C < C2.). 

• A large time region (usually t = 0(t^)) where C ^ I and X{C) < 1. In the same 
region the correlation function often satisfies an aging relation i.e C{t, t^) depends only 
on the ration s = t/t^ in the region where both t and are large: C(t, t^) ^ C^it/t^). 



In the simplest non trivial case, i.e. one step replica symmetry breaking 
function X{C) is piecewise constant, i.e. 

X{C) =m for C G /, X{C) = 1 for C ^ /. 



the 



(59) 



One step replica symmetry breaking for glasses has been conjectured in ref. [0. 

In all known cases in which one step replica symmetry holds, the quantity m vanishes 
linearly with the temperature at small temperatures. It often happens that m = 1 at T = Tc 
and m{T) is roughly linear in the whole temperature range. 
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7.2 Stochastic stability 



Stochastic stability is a property which is vahd in the mean field approximation; it is however 
possible to conjecture that is valid in general also for short range models. It has been 
introduced quite recently |3^, ^ |3^ and strong progresses have been done on the study 



of its consequences. 

In order to decide if a system with Hamiltonian H is stochastically stable, we have to 
consider the free energy of an auxiliary system having the following Hamiltonian: 

H + e^'^HR. (60) 

If the average (with respect to H^i) free energy is a differentiable function of e (and the 
limit volume going to infinity commutes with the derivative with respect to e), for a generic 
choice of the random perturbation inside a given class and e near to zero, the system is 
stochastically stable. 

The definition of stochastic stability may depend on the class of random perturbations 
we consider. Quite often it is convenient to chose as a random perturbation an infinite range 
Hamiltonian, e.g. 

Hr = Y^ Ji,k,icri<7k<Ji (61) 

i,k,l 

where sum runs over all the points of the system and the J's are random variables with 
variance 

In the nutshell stochastic stability tell us that the Hamiltonian H does not has any special 
features and that it properties are quite similar to those of similar random systems. 

Although it seems quite natural, stochastic stability has quite deep consequences. For 
example we could consider a system in which there are many equilibrium states, labeled 
by a, and the overlaps among two generic configurations of the same state and of two 
different states are respectively qm and qm, the free energies of the different states are 
uncorrelated. . . The situation would be quite similar to the one described by one step replica 
symmetry breaking. However we may not specify the form of the probability distribution 
of the free energies which is characterized by a function V{f) which a priori may have an 
arbitrary shape. 

It is a simple computation to verify that stochastic stability implies that the probability 
distribution of the free energies {V{f)) must have the form given in eq. (|38|) with an appro- 
priate choice of m. The most dramatic effect of stochastic stability is to link the behaviour 
of the function V{f) in the region of large / (where a large number of states do contribute) 
to the low / behaviour, which controls the distribution of the states which are dominant in 
the partition function. 

We have seen that stochastic stability strongly constraints the properties of the systems 
and many of the qualitative results of the replica approach can be derived as mere conse- 
quences of stochastic stability. Stochastic stability apparently does not imply ultrametricity, 
which seems to be an independent property |3^. This independence problem is still open as 
far has the only explicitly constructed probabilities distribution of the free energies of the 
states are ultrametric. 



7.3 Numerical simulations 

Let us consider the case of spin glasses at zero magnetic field (in this case the replica 
symmetry is fully broken |^). The natural variable to consider is a single spin {A = (jj)). 
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In this case the correlation C{t, t^) is equal to the overlap among two configurations at time 
t and t^: 

C{t,t^) = — . (62) 

The response function is just the magnetization in presence of an infinitesimal magnetic field. 
In this case the situation is quite good because there are reliable simulations for the system 
at equilibrium [^. 

In fig. (1) (taken from |22||) we plot the prediction for the function R versus C, obtained 



at equilibrium (i.e. using the equilibrium probability distribution of the overlaps, P{q)) by 
means of a simulation of a 16^ lattice using parallel tempering The simulation 

involves the study of 900 samples of a L = 16 lattice. 



During the off-equilibrium simulations p2[ in a first run without magnetic field the au- 



tocorrelation function has been computed. In a second second run from t = until t = 
the magnetic field is zero and then (for t > t^) there is an uniform magnetic field of small 
strength ho- The starting configurations were always chosen at random (i.e. the system is 
suddenly quenched from T = oo to the simulation temperature T). 

In fig. (1) there are the results of the off-equilibrium simulations p2| where = 10^ 
and tyj = 10^, with a maximum time of 5 ■ 10^ Monte Carlo sweeps. The lattice size in was 
64, and T = 0.7 (well inside the spin glass phase, the critical temperature is close to 1.0). 
We plot the response function R times T (in this case R is equal to m/ho) against C(t,t^). 
We have plotted also a straight line with slope —1 in order to control where the FDT is 
satisfied. Finally we have plotted two points, in the left of the figure, that are obtained with 
the infinite time extrapolation of the magnetization. 

The agreement among the absolute theoretical predictions (no free parameters) coming 
from the statics and the dynamical numerical data is quite remarkable. These data show 
the correctness of the identification of the functions x of the statics and X of the dynamics. 

Let us now go to the case of glass forming materials. I will present the data for binary 
mixture of soft spheres |^ . Theoretically there have been many speculations the glass tran- 



sition is described by one step replica symmetry breaking [|^, |S^, 0] . Here the equilibrium 
properties are no so well known as in spin glasses, although there are some evidence that 
the homologous of the function P{q) is non trivial PB[. On the other side, as we shall see. 



off-equilibrium simulations |^ show that the function X{C) seems to be given by the one 
step formula (|59|) with an approximate linear dependence of m on the temperature. 

We consider a mixture of soft particles of different sizes. Half of the particles are of type 
A, half of type B and the interaction among the particles is given by the Hamiltonian: 

+ . (63) 

where the radius (a) depends on the type of particles. This model has been carefully studied 
in the past [|l|, 0. The choice ctb/cta = 1-2 strongly inhibits crystallisation and the system 
goes into a glassy phase when it is cooled. Using the same conventions of the previous 
investigators we consider particles of average radius 1 at unit density. It is usual to introduce 
the quantity F = For quenching from T = oo the glass transition is known to happen 



around F^, = 1.45 41 



The best quantity we can measure to evidenziate off-equilibrium effects is the diffusion 
of the particles: 

Ht,tw) = ■ • (64) 
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The usual diffusion constant is given hy D = Yrnvt^^ A{t,t.u}) /t. 

The other quantity we measure is the response to a force. At time we add to the 
Hamiltonian the term e f ■ x^, where / is vector of squared length equal to d = 3 and we 
measure the response 



de 



£=0 



for sufficiently small e. The usual fluctuation theorem tells that at equilibrium (3A'^'^{t) = 

In the following we will look for the validity in the low temperature region of the gener- 
alized relation f]X{A) = dR/dA. This relation (with X 7^ 1) can be valid only in the region 
where the diffusion constant D is equal to zero. Strictly speaking also in the glassy region 
D 0, because diffusion may always happens by interchanging two nearby particles [D is 
different from zero also in a crystal); however if the times are not too large the value of D 
is so small in the glassy phase that this process may be neglected in a first approximation. 

The simulations we present are done using a Monte Carlo algorithm, which is a discretized 
form of a Langevin dynamics. In fig. ^we show R versus (3 A at t.^ = 2048 and = 8192 for 
r = 1.6 and t < At^ at N = 66. We also show the data for = 2048 at iV = 130. We do not 
observe any significant systematic shift in this plot among three data sets. We distinguish 
two linear regions with different slope as expected from one step replica symmetry breaking. 
The slope in the first region is compatible with 1, as expected from the FDT theorem, while 
the slope in the second region is near 0.62. Also the data at different temperatures for all 
values of r > 1.5 show a similar behaviour. The value of R, in the region where the FDT 
relation does not hold, can be very well fitted by a linear function of A as can be seen in fig. 

The region where a linear fit (with m < 1) is quite good corresponds to t/t^ > 0.2. 

The fitted value of m = dR/d{l3A) is displayed in fig. (3). When m becomes equal to 
1, the fluctuation-dissipation theorem holds in the whole region and this is what happens at 
higher temperatures. The straight line is the prediction of the approximation m(T) = T/T^, 
using Fc = 1.45. 

All the results are in very good agreement with the theoretical expectations based on 
our knowledge extracted from the mean field theory for generalized spin glass models. The 
approximation m(T) = T/Tc seems to work with an embarrassing precision. We can conclude 
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Figure 5: The quantity m = as = 2048 as function of the temperature. The straight 
hne is the prediction of the approximation m(T) = T /T^. 

that the ideas developed for generahzed spin glasses have a much wider range of application 
than the models from which they have been extracted. It likely that they reflect quite 
general properties of the phase space and therefore they can be applied in cases which are 
very different from the original ones. 

7.4 Some suggestions for planning an experiment 

The most interesting development would be to measure experimentally the function X both 
in spin glasses and in structural glasses. Clearly the most difficult task is the measurement 
of the fluctuations. In spin glasses it is clear how it should be done: the measurement of 
the thermal fluctuations of the magnetization is a delicate, but feasible experiment. In the 
case of structural glasses some ingenuity is needed in planning the experiments. (A open 
interesting possibility would to do the measurements in the case of rubber, where a transition 
with similar characteristics should take place.) 
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